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ABSTRACT 

We investigate the origin and fate of the recently discovered gas cloud G2 close to the Galactic Center. 
Our hydrodynamical simulations focussing on the dynamical evolution of the cloud in combination 
with currently available observations favor two scenarios: a Compact Cloud which started around the 
year 1995 and a Spherical Shell of gas, with an apocenter distance within the disk(s) of young stars 
and a radius of a few times the size of the Compact Cloud. The former is able to explain the detected 
signal of G2 in the position-velocity diagram of the Br7 emission of the year 2008.5 and 2011.5 data. 
The latter can account for both, G2's signal as well as the fainter extended tail-like structure G2t 
seen at larger distances from the black hole and smaller velocities. In contrast, gas stripped from 
a compact cloud by hydrodynamical interactions is not able to explain the location of the detected 
G2t emission in the observed position-velocity diagrams. This favors the Spherical Shell Scenario and 
might be a severe problem for the Compact Cloud as well as the so-called Compact Source Scenario. 
From these first idealized simulations we expect a roughly constant feeding of the supcrmassive black 
hole through a nozzle-like structure over a long period, starting shortly after the closest approach in 
2013.51 for the Compact Cloud. If the matter accretes in the hot accretion mode, we do not expect 
a significant boost of the current activity of Sgr A* for the Compact Cloud model, but a boost of the 
average infrared and X-ray luminosity by roughly a factor of 80 for the Spherical Shell scenario with 
order of magnitude variations on a timescale of a few months. Assuming that a part of the gas is 
accreted in cold disk mode, even higher boost factors can be reached. The near-future evolution of the 
cloud will be a sensitive probe of the conditions of the gas distribution in the milli-parsec environment 
of the massive black hole in the Galactic Center. 
Keywords: accretion - black hole physics - ISM: clouds - Galaxy: center 



1. INTRODUCTION 

The Galactic Center harbors a nuclear supermassive 
black hole (SM BH). Due to its proximity, its mass of 
4.31 X 10^ Mp, jGillessen et af 



2009 



Ghez et al. 



2008 



Genzel et al. 2010) has been derived from direct long- 



Montero-Castano et al. 2009) and a shell-like, elliptical 
non-thermal structure with a size of roughly 6x8 pc 
(called Sgr A E ast, Ekers fc Lynde n-BcU 197ll [ Dowrie s 
& Martin 119711 lEkers et al. 1975; Macda et al. 2002^ 



term (w20 years) monitoring of stellar orbits in the sub- 
parsec-scale vicinity of the central compact nonthermal 
radio source Sgr A* - thought to be the radiative coun- 
terpart of the central SMBH. The uncertainties in the 
mass-estimate of the order of 10% are mostly due to the 
distance measurement. The parsec-scale environment 
of th e center i s compr ised of seven main c onstituents 
(e. 



I'br a more detailed review of the Galactic Center envi 
ronment we refer to 



Recently, a gas c 



Genzel etal.| ( |2010 ). 
oud has been detected on its way 

It 



2012) 



Melia & Falcke 2001 Genzel et al.|2010|): an atmo- 



^ sphere ot d ittuse hot gas ( e. g . Bagaiioff et al.|2003 Yuan 



et al.|[2003 se e also Sect. 3.1), a cluster ot evolved stars 



( [Genzel et al. 



stars ranging 



20031, one or two warped disks of young 



jetween 0.04 pc and 0.5 pc (Genzel et al. 



20031 [Levi n fc Beloborodwl [2003| [Bartko e t al.[[:^0i01 



Lu et al. '2009; Alig et al. 20TTt, the B-stars in the 



central half-parsec ( ,Gillessen et al. 2009, ), spiral shaped 
thermal gas streanis in the inner 2 pc, which reaches 
down to 0.17 pc (called Sgr A West, " Ekers fc Ly nden-Bell' 
1971; Downes & Martin 1971 ; Ekcrs et al.| |iij83; Yusef- 



TaSehTk Wardle|[1993[|ZEao~et al.[[2009D, a molecular 



dusty ring ranging trom 2pc to 5 pc ( Genzel et al.,,1985[ 
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towards the Galactic Center (Gillessen et al 
consists of a two-component system: (i) the high surface 
brightness compact component G2 and (ii) a lower sur- 
face brightness cone-like structure, which we will refer 
to as G2t fro m now on, as it was o riginally described as 
the tail of G2 ( Gillessen et al.|2(3T2 ) . The whole complex 
(G2 plus G2t) will be called eG2, the extended C2 cloud 
complex, in o rder to avo id confusion. Traced back on 
archival data, Gillessen et al.[ ( [2012 ) were able to deter- 
mine an accurate orbit tor G2 (I'ig. ll) , which has a very 
high eccentricity (e = 0.9384). It will reach its closest 
approach of only 4 x 10^^ cm to the SMBH in 2013.51 
(corresponding to 3100 Schwarzschild radii). Clear evi- 
dence for ongoing tidal velocity shearing and stretching 
of G2 within the last couple of years has been detected, 
confirming its extended nature. This provides an inter- 
esting challenge for hydrodynamical simulations that aim 
at reproducing the currently visible tidal shearing and 
that can be used to predict the future tidal evolution 
over the next couple of years. G2 comprises of a multi- 
phase medium with the bulk of its mass being ionized 
gas {Mq2 ~ 1.7 X 10^^ /vg, with /y being G2's volume 
filling factor, which we assume to be unity throughout 
this paper). Ionization equilibrium with the radiation 
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Figure 1. Orbits of our two scenarios: Tlie black line represents tlie orbit of G2, as determined in'Gillessen et aL](|2012[ l, representative 
for the Compact Cloud Scenario and the red line denotes the orbit of eG2 as determined in this work, representative tor the Spherical Shell 
Scenario. The location of the central black hole is denoted by the plus sign and the dashed line denotes the inner rim of the disk(s) of 
young stars. The sizes of the shell and the cloud are shown by the thin black lines. 

field of tlie nearby massive stars is expected to keep the 
temperature of this component at around 10'* K. Addi- 
tionally, it contains a small amount of warm dust with a 
temperature of ~ 550 K. Only very little is known about 
G2t, but its mass could be substantially larger than the 
mass estimate for G2 itself. 

Several questions are of immediate theoretical interest: 

• Which physical process formed the cloud? 

• What determines the physical cloud properties 
(mass, size, thermodynamic state) and the dynam- 
ics? 

• Can it be used to probe the direct vicinity of 
SgrA*? 

• What is the fate of the cloud in the next years? 

• How does it influence the activity of Sgr A*? 

In IBurkert et al. (2012) (paper I) we mainly tackled 
the first three of these questions mostly with the help of 
analytical considerations. We especially concentrated on 
possible formation scenarios and evolutionary paths of 
G2. Two basic models proved to be the most successful: 

• The Compact Cloud Scenario (CC) 

A cloud with uniform density structure can explain 
the evolution of the observed position-velocity 
(PV) diagrams of G2. 

• The Spherical Shell Scenario (SS) 

A tidally disrupted, initially spherical shell of gas 
with a much larger radius and mass than G2 can 
account for both, the G2 and G2t component in 
the position-velocity diagram. 

In paper I we also found that the cloud is a sensitive 
probe of the hot accretion zone of the SMBH and the cur- 
rently available observations of G2 are in good agreement 
with models of captured, shock-heated stellar winds. We 
were able to derive constraints on the mass of G2, which 
is above the critical mass below which clouds would evap- 
orate before reaching the presently observed radius, and 
which is small enough for the sound crossing timescale 



at apocenter to be of order its orbital period. Analytical 
estimates of various other time scales showed that many 
physical processes might contribute to the evolution of 
the cloud. Therefore, detailed numerical simulations are 
needed, which is the topic of this publication. The wealth 
of available observational data (now and in future) as well 
as its proximity make this event an ideal tool to study 
gas accretion physics in the vicinity of SMBHs. It will 
also give us very valuable insight into the evolution and 
the activation of galactic nuclei in general. 

An alternative model to our Compact or Spherical 
Cloud Scenario - called th e Compact Source Sceriario in 
paper I - was proposed by [Mu rray-Clay fc Loeb| (20111 
(see also Miralda-Escude 2012): a dense, proto-planetary 
disk, which is bound to its parent low-mass star. The star 
is assumed to be scattered roughly 100 years ago from the 
observed disk(s) of young stars. Far away from SgrA*, 
the protostellar gas disk remains bound to the star and 
resists tidal stretching. When approaching the SMBH, 
it looses mass through photo-evaporation due to the ex- 
treme flux of ionizing and FUV ph otons in the Galactic 
Cente r and due to tidal stretching. [Murray-Clay fc Loeb| 
(20111 analytically estimate that tnis results in a cloud 
with similar properties as observed. 

The main idea of this paper is to present a first set 
of idealized hydrodynamical simulations and to test the 
influence of the hot atmosphere on the late-time evo- 
lution of such a cloud for our two best-fit cloud infall 
models as already described in paper I, namely the Com- 
pact Cloud Scenario and the Spherical Shell Scenario. 
We discuss simple test particle simulations in Sect. [2] be- 
fore we describe the model we employ for the Galactic 
Center environment in Sect. |3| and review some of the 
complicated physical processes in Sect.[4| The results of 
our hydrodynamical simulations are presented in Sect.[5J 
followed by a discussion (Sect.[6| and finally the conclu- 
sions (Sect.[7]). We assume a Cfalactic Center d istance of 
8.33 kpc, w here 1" = 40mpc = 1.25 x 10^^ cm ( [Gillessen 
e raL][2009| . 



2. TEST PARTICLE SIMULATIONS 

In order to find initial conditions for our hydrodynam- 
ical studies, which are able to reproduce current obser- 
vations, we calculate a grid of test particle simulations in 
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two dimensions, varying a large number of the parame- 
ters, as described below. This is d one for both models of 
origin as discussed in Sect. 2.1 and 2.2 Such kind of sim- 
ulations are very fast and hence a large number of test 
calculations can be done on a short timescale (approxi- 
mately 16 000 simulations for this paper), allowing for a 
huge scan of the involved parameters. The test particle 
simulations are ru n with the N-body code REBOUND 
( Rein fc Liu||2012 ) and 100 test particles are distributed 
along the boundary of the cloud or ring, restricted to 
the orbital plane. The particles are then evolved in the 
gravitational potential of the black hole and describe the 
evolution of the shape of the cloud when tidal forces com- 
pletely dominate. The best fit model is found by com- 
paring the simulated to the observed position velocity 
diagrams. We use these N-body simulations to roughly 
constrain the formation time and initial conditions of G2 
and eG2, respectively. The best fitting initial conditions 
are then used to study the detailed hydrodynamical evo- 
lution (Sect. [5]). 

2.1. Compact Cloud Scenario 

For the compact cloud scenario, we assume that the 
cloud started in situ in pressure equilibrium with the 
surrounding hot atmosphere on the observed orbit of G2 
and that G2's mass is knowrF] Then, the initial radius 
and density of the cloud are nxed along the orbit and we 
only vary the date of birth of the cloud. As estimated 
in paper I, hydrodynamical effects - especially ram pres- 
sure - play a minor role in the early evolution of the 
cloud in the compact cloud scenario. This means that 
we can directly compare the resulting distribution in the 
position-velocity diagram of the year 2008.5 and 2011.5 
with the observations and can expect to get a reasonably 
well determined initial condition. The goodness of fit is 
judged from the position and spread of the test particles 
in the PV diagram. We compare the minimum and max- 
imum values of the projected distance {(P^.^^^) as well 
as the line of sight velocity (w^™ ^^x) o^ the simulated 



distribution with the respective extrema of the observed 



1,000 








^ 








- 


0,100 


\ 




^ 


\ 


X 


^> 


\ / 


4 


0,010 






X 


M 


X 


\ .11, 


^ 


/, 


0,001 






10 


Ho 




m 


Hfl 0.5 Ho 


■ 



1940 1950 1960 1970 1980 1990 2000 

formotion time 

Figure 2. Possible formation times of the compact cloud models 
as inferred from test particle simulations for various densities of 
the atmosphere. The black line corresponds to our standard atmo- 
sphere, whereas the colored lines refer to multiples of the density 
as indicated in the plot. The rp value quantifies the goodness of fit 
with the obse rved position-velocity diagrams in 2008.5 and 2011.5 
(see Sect.|2.1|l. 

** The observationally determined mass actually depends on the 
unknown volume filling factor, which we assume to be one here. 



PV diagram [v. 



obs 



Jobs 



min,max' min,max 



) . For the latter we chose 



a representative contour line (0.8 times the maximum 
value of the observed PV diagram, shown by the two plus 
signs in the left two panels of Fig. 10 1. The quadratic 



sum of these relative deviations of th^008.5 and 2011.5 
observations with the respective simulated snapshots are 
summarized in an 77^-value, used to judge the goodness 
of the fit: 
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mm, max 
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(1) 



Fig. [2] displays the result of t his st udy for our assumed 
standard atmosphere (see Sect. 3.1). It shows a clear rf 
minimum around the year 1995.5. A study of different 
density normalizations of the atmosphere (Fig. l2| shows 
that the preferred origin of the cloud depends sensitively 
on the density distribution of the atmosphere. Increas- 
ing the atmospheric density by a factor of two shifts the 
starting time to the year 1988, whereas a factor of ten 
would necessitate a starting time of 1967 and a decrease 
of the density distribution by a factor of two yields a 
starting date of roughly 2000. This can be understood as 
the cloud sizes depend on the atmospheric pressure, when 
we assume pressure equilibrium between the cloud and 
the atmosphere. However, the values also depend on the 
choice of the limiting contour for the fitting procedure. 
Therefore, these studies only serve as rough estimates for 
possible origin dates of the cloud. To determine and de- 
scribe the future hydrody nam ical evolution, we will start 
the cloud in 1995.5 (Sect.[5l|). 

2.2. Spherical Shell Scenario 

As shown in paper I, the observations can also be inter- 
preted in such a way that G2 and G2t are parts of a larger 
cloud complex (eG2), which is already partly disrupted 
by tidal forces in the 2008.5 and 2011.5 observations. As 
these shells are typically a factor of five to ten larger 
than the compact clouds mentioned above, ram pressure 
effects cannot be neglected anymore. In order to find 
reasonable starting parameters for this model as well, 
we take the ram pressure effects roughly into account by 
offsetting the comparison values from the observed con- 
tour (compare to blue plus signs in panels c,d of Fig. 10 ). 
This offset is determined by comparing a number of hy- 
drodynamical simulations with test particle simulations. 
Thought to be the result of emission of stellar material, 
we assume that this shell started close to the apocenter 
of its elliptical orbit, which lies well within the disk(s) of 
young stars (see Fig.Il]), making it the most likely origin 
of the cloud. For this model, the total mass of the cloud 
is not fixed, as the mass of G2t cannot be constrained 
observationally very well at present. In the subsequent 
hydrodynamical modelling we then choose the density of 
the shell to reach approximate pressure equilibrium with 
the surrounding atmosphere. Again, we are only inter- 
ested in a rough match of the hydrodynamical realization 
with the observed PV-diagram in order to investigate the 
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Figure 3. Distribution of r^'^-values (see Sect. |2.l| of the spherical shell scenario as inferred from test particle simulations in the initial 
x-position vs. initial y-velocity plane. We assume the ring starts at arbitrary time at the apocenter, which is located on the negative y-axis. 
The three different panels refer to three exemplary initial ring radii: 1.00 X 10^^ cm (a), 1.18 X 10^^ cm (b) and 1.36 X 10^® cm (c). The rj^ 
value quantifies the proximity of the match with the observed position-velocity diagrams in 2008.5 and 2011.5 (see text for details of the 
fitting procedure). 

basic evolution of such a model. In this case, every time 
step is compared to the 2008.5 PV diagram and the stage 
three years later with the 2011.5 observation in order to 
find the best matching timestep, as the date of birth is 
not fixed in this model. The result of this fitting proce- 
dure with a grid of 25'^ test particle simulations is sum- 
marized in Fig. pi Shown is the distribution of r]^ values 
determined as described in Sect. 2J_ and equation [T] in 
a plane spun by the initial position on the negative x- 
axis and the initial velocity in y-direction. The panels 
display the distribution for three examples of the initial 



(Bower et al. 2003). In the idealized simulations pre- 
sented in this paper, we neglect the rotation and the ac- 
tual dynamics of the accretion flow and rescale the nor- 
malization of the temperature distribution accordingly, 
in order to yield hydrostatic equilibrium. This yields the 
following number density distribution riat and tempera- 
ture distribution Tat of the hot atmosphere: 



Tiat = 930 cm •^ /i 



/hot 



outer ring radius: 1.00 x 10 cm (panel a), 1.36 x 10 cm 
(panel c) and our best fit value of 1.18 x lO^^cm (panel b). 
The minimum ij^ valley is quite shallow and we take the 
global minimum of our simulated grid, which is a model 
starting at apocenter in the year 1927.2 at a distance 
from the black hole of 1.58 x 10^^ cm and an initial or- 
bital velocity of 125kms~^. 

3. MODELLING THE GALACTIC CENTER ENVIRONMENT 
3.1. The hot atmosphere 
X-ray observations of hot gas emission near SgrA'' 



Tat, = 1.2 X 10«K 



1.4 X lO^i^s 



(2) 
(3) 



( Baganoff et al. 



2003) are attributed to the atmosphere 



close to the HJVIBH, which can be explained as gas being 
shock-heated by the interaction of stron g stellar winds 
of th e surrounding massive stars (e. g. iKrabbe et al. 



1991 



2004 



|Melia| [19921 [Baganoff et~ar] [2003lnJVluno et a 



The observed very low bolometric luminosity of 
the order of 3 x 10~^ LEdd to gether with a Bon di ac- 
cretion rate of 10'^ M^yv-^ ( [Yuan et a f 2003 ) disfa- 
vor s tandard accretion disk theory ( bhakura fe Sunyaev] 
19731. They are rather consistent with so-called Hadia- 
tively Inefficient Accretion Flows (RIAFs), in which only 
a small fr action of the gravitatio nal energy is radiated. 
Following [Gillessen et al.[ (|2012), we adopt t he ADAF 
( Advection-Dominated ~Kccretio n Flow, e. g. Ichimaru 
19771 [Abramowi czet ar 'l995: Na rayan fc Yill994p model 

fOU3), which is an analytical approxima- 
: 



where r is the distance to the black hole, Ra, is its 
Schwarzschild radius, a « /3 ~ 1 and /hot ~ 1 is a fac- 
tor taking the uncertainty of the model into account. A 
mean molecular weight of /i = 0.6139 has been assumed, 
typical for a gas with solar metallicity. 

However, as the entropy S = T/ p^/^ decreases out- 
wards, this is not a stable equilibrium, but unstable 
to convection (see discussion in paper I), which we con- 
firmed with the help of test simulations: small pertur- 
bations start growing at the inner boundary surround- 
ing the black hole, which quickly leads to the forma- 
tion of unphysically high pressure bubbles. Due to the 
unstable nature of the atmosphere we do expect an in- 
homogeneous density distribution. However, setting up 
a real istic RIAF solution is a complicated task in itself 



e. g. Stone et al. 1999, Stone fc Pringle 2001[and Igu- 



menshcEev et al.,2003J and beyond the scope oi this hrst 



oF YuaiTet al 
tion tor such a 



idealized study presented in this article as it would re- 
quire the inclusion of many m ore physical processes, like 
e. g. magnetic fields (see Sect. 6.3 for a discussion of pos- 
sible effects), thermal conduction, etc. Some of these 
effects are already taken into acco unt (not com pletely 
" ' ' ' (2003). 



self-consistently) in the models of Yuan et al. 
We are rather interested in the interaction ot the"~cioud 
with a smooth atmosphere, representing an average dis- 
tribution in concordance with available observations. To 
rotation measure data derived from radio observations this end, we artificially stabilize it and keep the atmo- 



?I70^ 



solution , in concordance with avail- 
able Chandra observations (Baganoff et al. 2003) and 
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sphere's density and temperature distribution fixed with 
time. This is done by additionally evolving a passive 
tracer field (0 < ir < 1), which allows us to distinguish 
between those parts of the atmosphere whic h ha ve in- 
teracted with the cloud {tr > 10~^, see Sect. 6.2) from 



those which changed due to the atmosphere's mhcrent 
instability {tr < 10"'*). Those cells fulfilhng the latter 
criterion are reset to the values expected in hydrostatic 
equilibrium. We discuss numerical is sues and proof that 
this is a sensible assumption in Sect. |6.2| below. 



3.2. The initial conditions of the cloud models and the 
numerical treatment 

We use the two best-fit models presented in Sect. [2] as 
basis for our hydrodynamical exploration of the origin 
and fate of the cloud complex eG2. In these first simula- 
tions, we generally set the adiabatic index F = 1, but also 
present the comparison to a purely adiabatic simulation. 
This seems to be a reasonable assumption, as the tem- 
perature structure of the atmosphere is supposed to be 
given by adiabatic heating of th e accretion flow itself and 
stellar processes (see Sect. 3.11 and the temperature of 
the cloud material is expected to be set by photoioniza- 
tion e quilibrium in the radia tion field of the surrounding 
stars (Gillessen et al. 2012) during the early evolution 
of the cloud. In order to explain the observed position- 
velocity diagram, two basic scenarios have been identified 
in paper I. Both are based on the simple assumption that 
the cloud started with spherical shape and in pressure 
equilibrium with the atmosphere on an orbit compatible 
with the observed positions on the sky. In the first set 
of mo dels we follow the interpretation of |Gillessen et ah 
(2012), who discern the observations into cloud head and 
tail emission. Our second set of models interprets the 
observations as being the result of the disruption of a 
larger cloud, which - for concordance with observations 
- necessitates the cloud to be a spherical shell. 

Our two-dimensional simulations are done within the 
orbital plane of the cloud, with the major axis parallel 
to the x-axis and the apocenter of the orbit on the nega- 
tive x-axis of our fixed Cartesian coordinate system with 
a spatial resolution of 7 x 10^'^cm. The computational 
grid of our standard model ranges from —1.3 x lO^^cm 
to 1.2 X lO^^cm in x-direction and —6.2 x lO^^cm to 
2.5 X lO^^cm in y-direction. The cloud is on a clock- 
wise orbit with the major axis parallel to the x-axis and 
the pericenter of the orbit on the positive x-axis. The 
black hole is located at the origin of our coordinate sys- 
tem. As discussed above we use a passive tracer field in 
order to trace the evolution of the cloud (and to stabilize 
the atmosphere) and initially assign the cloud a value of 
1 and the atmosphere and let it evolve. 

The hydrodynamical equations are integra ted with 



PLUTO, version 3.1.1 ( [Mignone etaT 
simulations shown in this article 



2007D . For aU 
we use the two-shock 



Riemann solver, do a parabolic interpolation and employ 
the third order Runge-Kutta time integration scheme. A 
cartesian two-dimensional grid is used and the boundary 
conditions in all directions are set to the values expected 
for hydrostatic equilibrium of the hot atmosphere, en- 
abling outflow but no inflow. We introduce an additional 
spherical boundary surrounding the central black hole. 
All gas flowing into this part of the integration domain 
is removed and counted as being accreted. All values in- 



side this boundary are set to the atmospherical values, 
in order to mimick the pressure exerted by this region, 
which is necessitated by the steep pressure gradient. 

We neglect thermal conduction, magnetic flelds as well 
as feedback from the central source for the sake of sim- 
plicity. A summary of the used and varied parameters is 
given in Table[T] and Table[2] respectively. 

4. CLOUD PHYSICS 

The components of the Galactic Center environment 
which have the largest impact on the gas cloud are the 
gravity of the central SMBH and the hydrodynamical in- 
teraction with the hot atmosphere. Therefore, we restrict 
our discussion to the physical effects related to these 
two components. The relevant physical processes of the 
cloud/atmosphere system are discussed in great detail in 
paper I. Here, we only briefly summarize those important 
for the interpretation of our hydrodynamical simulations 
and estimate the timescales of the various processes for 
different phases of the evolution of the cloud. 

4.1. Tidal deformation 

Tidal forces lead to a stretching in direction of the 
black hole and a squeezing perpendicular to it. As G2's 
orbit has a very high eccentricity, the cloud gets dis- 
tributed along the orbital path. The acceleration of this 
stretching is given by 



|atf| 



2ArG 



M, 



BH 



(4) 



where Ar refers to the width of the cloud in direction 
of the black hole, G is the gravitational constant, A/bh 
is the mass and r the distance to the SMBH. Therefore, 
a change of the cloud elongation of the order the width 
of the cloud by tidal forces is expected on a timescale of: 



12 Ar 



r3/2 



Ttf 



|att| VGM^ 



10.60 yr 



JBH 



1016 cm 



3/2 



(5) 



which varies between 27.6 yr at the 1995.5 position of 
the cloud and 0.5 yr at the pericenter. 

4.2. Ram pressure compression and deceleration 



Table 1 

Observed and calculated parameters of the cloud G2 and the 
Galactic Center region. 



Parameter 



Value 



mass of the SMBH 

Galactic Center distance 

mass of the cloud G2'' 

pericenter distance 

time of closest approach 

orbit eccentricity 

semi-major axis 

orbital period 

inclination of ascending node 

position angle of ascending node 

longitude of pericenter 



A/bh 


4.31 X 10"^ Mq 


i?.o 


8.33 kpc 


JTlcloud 


1.7 X 102**g 


^pcri 


4 X 10^5 cm 


^pcri 


2013.51 


e 


0.9384 




6.49 X IQiScm 


^orbit 


137.77 


I 


106.55 


n 


101.5 


LJ 


109.59 



Note. — Parameters of G2 and the Galactic Center obtained or 
derived from observations in Gillessen et al. (2012). (a) Through- 
out the paper we assume a volume hlling factor of one. 
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Table 2 

Parameters of the (hydrodynamical) simulations. 





to" 


Pcloud "^ 


Kcloud'^ 




J/ini" 


V? .^ 


v".^ 


rh 


/hot' 


AxJ 


tr^ 




yr AD 


10-"gcm-3 


10^5 cm 


lO^^cm 


10l«cm 


kms ^ 


kms ' 






10^3 cm 




CCOl 


1995.5 


6.21 


1.87 


-7.22 


2.21 


794.59 


48.45 




1 


7.0 


10-4 


CC02 


1944.6 


2.24 


2.63 


-12.59 


0.0 


0.0 


167.29 




1 


7.0 


10-4 


CC03 


1944.6 


223.64 


0.57 


-12.59 


0.0 


0.0 


167.29 




1 


7.0 


10-4 


CC04 


1995.5 


6.21 


1.87 


-7.22 


2.21 


794.59 


48.45 


5/3 


1 


7.0 


10-4 


CC05 


2000.0 


4.33 


2.11 


-5.99 


2.24 


953.00 


-5.41 




0.5 


7.0 


10-4 


CC06 


1988.0 


8.57 


1.68 


-8.87 


2.03 


605.85 


98.85 




2.0 


7.0 


10-4 


CC07 


1995.5 


6.21 


1.87 


-7.22 


2.21 


794.59 


48.45 




1 


14.0 


10-4 


CC08 


1995.5 


6.21 


1.87 


-7.22 


2.21 


794.59 


48.45 




1 


3.5 


10-4 


CC09 


1995.5 


6.21 


1.87 


-7.22 


2.21 


794.59 


48.45 




1 


7.0 


lo-'' 


CCIO 


1995.5 


6.21 


1.87 


-7.22 


2.21 


794.59 


48.45 




1 


7.0 


10-2 


SSOl 


1927.2 


1.42 


11.80' 


-15.80 


0.0 


0.0 


125.00 




1 


7.0 


10-4 



Note. — CC refers to simulations of the Compact Cloud scenario and SS to those of the Spherical Shell scenario. Simulations CC02 
and CC03 are only discussed in paper I. (a) Start time of the simulation, (b) initial density of the cloud, (c) initial radius of the cloud, 
(d) initial x-position of the cloud, (e) initial y-position of the cloud, (f) initial x-velocity of the cloud, (g) initial y-velocity of the cloud, 
(h) adiabatic index, (i) uncertainty factor of the density distribution of the hot atmosphere, (j) cell size in x- and y-direction, (k) tracer 
threshold value, (1) the shell has a thickness of 25% of this outer radius. 



Ram pressure of the surrounding medium has two ef- 
fects: (i) LocaUy, it changes the cloud structure by com- 
pressing the upstream part of the cloud and (ii) glob- 
ally it decelerates the cloud, leads to angular momentum 
exchange with the atmosphere and finally accretion of 
matter towards the center. In paper I we analytically es- 
timate the loss of kinetic energy for our compact cloud 
model to be of the order of a few percent under the 
assumption of a cloud evolving in pressure equilibrium. 
In realistic hydrodynamical simulations, the situation is 
much more complex as we will show in detail in Sect. [5] 

4.3. Thermal conduction and evaporation 

Thermal conduction l eads to the evaporation of parts 
of the cloud. Following |Cowie fc McKee" (1977) we find 
an evaporation time scale of 



Tcvap = 81yr 



10 r. 



15/8 



Re 



10 



15 



13/8 



cm 



1.7 X 1028 g 



(6) 



in the so-called saturated limit for a spherical cold gas 
cloud immersed in the hot atmosphere described in equa- 
tion[2]and[3| Assuming that the cloud stays spherical and 
in pressure equilibrium along the whole orbit, the evap- 
oration time scale changes between 97 yr for the cloud 
at its initial distance in the year 1995.5 and 10 yr when 
the cloud reaches the pericenter. However, this neglects 
tidal effects and hydrodynamical instabilities, which fi- 
nally lead to the disruption of the cloud. Thereby, the 
cloud increases in size and the evaporation time scale 
might decrease. Hence, evaporation effects get impor- 
tant mostly during the late evolution of the cloud. 

As was shown in paper I and briefly discussed here, all 
of the mentioned physical processes are acting on simi- 
lar timescales close to the orbital timescale of the cloud. 
This is especially the case for the late-time evolution and 
necessitates the use of numerical hydrodynamical simu- 
lations to determine the fate of the cloud. 



5. RESULTS OF THE HYDRODYNAMICAL SIMULATIONS 

5.1. The evolution of the density distribution in the 
Compact Cloud Scenario (CCOl) 

The temporal evolution of the two-dimensional gas 
density distribution of model CCOl is shown in Fig. [4J 
Overlayed are the results of a simulation in which we 
placed test particles at the boundary of the initial cloud 
setup and let them evolve under the influence of the grav- 
itational potential alone. They serve as an indicator for 
the differences in the evolution due to hydrodynamical 
effects. We start with a phenomcnological description of 
the cloud evolution and back our findings in the course of 
the paper with several analysis tools. The initially spher- 
ical cloud begins its journey at the year 1995.5 position 
on G2's observed best-fit orbit. The early phase towards 
the pericenter which is reached in 2013.5 is dominated by 
tidal forces (panels a to e). Due to the highly eccentric 
orbit, the cloud gets stretched along the orbital direction. 
At the same time, the surrounding density and the in- 
crease in orbital velocity lead to an increase of ram pres- 
sure forces, resulting in a compression of the cloud head 
which can be inferred from the deviation of the density 
distribution from the test particle model (panel d) . The 
shear flow at the boundary of the cloud leads to stripping 
of gas and the formation of a long turbulent wake, which 
develops typical rolls of the Kelvin- Helmholtz instability 
(KHI, panel e). At the time of the pericenter passage, 
the cloud has developed into a long spaghetti-like fila- 
ment more or less aligned with the test particle orbit of 
the cloud. This elongation of the cloud in combination 
with the large velocity and density difference (along G2's 
orbit) between the cloud and the atmosphere causes a 
rapid growth of the KHI, well visible by the meandering 
shape of the cloud (panel f). In the subsequent evolu- 
tionary phase, a number of effects lead to the buckling of 
the cloud in direction of orbital motion (panel g): First, 
the pure gravitational effect, leading to a deceleration of 
the front while the back still has a high velocity. This is 
the dominant effect, as the horizontal width of the test 
particle distribution almost coincides with the horizontal 
width of the gas cloud. Second, this is the region of the 
orbit, where the atmosphere is the densest. Hence, strong 
ram pressure forces shape the head of the cloud and lead 
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Figure 4. Density evolution of the compact cloud model CCOl. Overlayed as dotted white lines are the positions of test particles initially 
located at the cloud boundary. The axis labels are given in mpc (milli-parsec). Mind the different scalings of the panels. 



to additional compression along the orbit and elongation 
perpendicular to the orbital motion. Increasing the cross 
section, ram pressure forces get more and more impor- 
tant. A boundary layer builds up at the front of the 
cloud, which looses angular momentum and cloud mate- 
rial is able to accrete towards the center, forming a num- 
ber of filaments pointing towards Sgr A* (panel g) . Their 
number and location is dictated by the earlier growth of 
the Kelvin-Hclmholtz instability mainly during the peri- 
center passage. In the course of the following evolution, 
the rear filaments - being in the shadow of the leading 
one - merge with the one in front and form a nozzle- 
like structure of matter streaming towar ds Sg r A* and 

panel h) . 
by hydro- 



5.3 



causing an increased accretion rate (Sect. 
The later evolution of the cloud is dominated 
dynamical effects in our simulations: The large effective 
area of the cloud makes it susceptible to hydrodynami- 
cal instabilities and ram pressure. The phase during and 
shortly after the pericenter passage is also the time when 
for the first time, the cloud is able to mix with the atmo- 
sphere. Due to the large temperature difference between 



cloud and atmosphere and the moderate density contrast 
in this region, this mixing is able to increase the cloud 
temperature to approximately 10^ — 10® K. At the end of 
this mixing process, about 1% of the cloud is made up of 
atmospheric gas. 

Whereas hydrodynamical instabilities shape the outer 
region of the cloud, ram pressure interaction with the 
atmosphere enables further deceleration and accretion of 
gas from the boundary layer. However, it is very dif- 
ficult to make detailed predictions of the late evolution 
of the gas with the help of our idealized simulations, as 
by then, several other effects become important. Ther- 
mal conduction will be able to act on the much increased 
surface area of the cloud and will evaporate it outside- 
in and the detailed thermodynamical treatment might 
as well significantly change the morphology and accre- 
ti on b ehavior in this evolutionary phase (see discussion 
in 



5.5). 
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Figure 5. Density evolution of tiie spherical shell model SSOl. Overlayed as dotted white lines are the positions of test particles initially 
located at the outer ring boundary. The axis labels are given in mpc (milli-parsec). Mind the different scalings of the panels. 



5.2. The evolution of the density distribution in the 
Spherical Shell Scenario (SSOl) 

Density snapshots showing the evolution of the spher- 
ical shell scenario (simulation SSOl) are displayed in 
Fig. [5) The basic physical processes acting on the orbit 
untilthe pericenter passage are very similar to the evo- 
lution of the compact cloud. Firstly, the ring squeezes 
into a drop-shaped morphology due to the interaction 
in the gravitational potential of the black hole (panel a) . 
Secondly, the hydrodynamical interaction with the hot 
atmosphere leads to the formation of a turbulent wake 
behind the shell and at the ring border, KHI start to grow 
(panel b). In a three-dimensional simulation, the struc- 
ture would correspond to a hollow conical morphology, 
with one major focal point at the very front of the cloud. 
This is in good agreement with what is seen in obser- 
vations: The focal point then corresponds to the actual 
cloud G2 and the conical structure is the diffuse tail G2t, 
which indeed shows signs of a conical shape. In reality, 
the conical appearance is irregular and not continuous 
as the spherical shell is expected to have an inhomo- 
geneous density distribution initially. At pericenter, the 
front part of the ring - which we interpret as the cloud G2 
- has already collapsed into a thin, spaghetti-like s hape . 
In contrast to our compact cloud model (see Sect. 5.1) 



the thin filament shows already the typical disturbances 
expected for the KHI and hence a larger effective cross 
section compared to the compact cloud (panel b). Be- 
ing very close to the center of the hot atmosphere, the 
steep density gradient causes a rapid increase of the ram 
pressure, decelerating the front part of the ring. The 
KHI disturbances already present before pericenter pas- 
sage grow on short timescale and finally lead to a fast 
disruption of the ring (panel c) . It efficiently exchanges 
angular momentum with the dense inner atmosphere and 
a filamentary disk- like structure is able to form (panel d) . 
This leads to a much higher mass accretion rate through 
the inner boundary (see Sect. 5.3 ). We refer to Sect.lHlfor 



a critical discussion of the infiuence of numerical effects 
as well as our idealized treatment of the atmosphere. 

5.3. Angular momentum redistribution and mass 
accretion towards the center 

The temporal evolution of the binding energy of the 
cloud (selected with the usual tracer threshold limit) is 
shown in Fig. |6] for the two models CCOl and SSOl. As 
already discussed in paper I, the loss of kinetic energy due 
to ram-pressure in the early evolutionary phase is minor 
and the binding energy stays more or less constant until 
the first part of the cloud reaches the dense inner part 
of the hot atmosphere around the year 2012/2013, where 
we see a strong increase of ram pressure deceleration. As 
a consequence, the cloud gets bound more strongly to 
the SMBH. Ram pressure forces increase with the effec- 
tive cross section of the cloud, they rise sharply shortly 
after the pericenter passage of the CC model, when the 
cloud compresses in direction of motion and expands per- 




2000 2010 2020 2030 2040 2050 2060 
time [yr] 

Figure 6. Binding energy of the cloud/ring in units of the initial 
binding energy. 
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pendicular to the orbital direction. A small fraction of 
the gas in the mixing zone of the cloud and atmosphere 
then carries angular momentum outwards, enabling part 
of the cloud to be accreted towards the center. When 
this outward moving gas has become very diluted and 
mixed with the atmosphere, the angular momentum is 
removed from the simulation due to our treatment to 
stabilize the unstable atmosphere. After the pericenter 
passage a dilute high pressure bubble forms (see Fig. |4 
panel i), which rises buoyantly, thereby pushing some of 
the gas outward and the remaining gas gets less bound to 
the SMBH (Fig.|6f. Model SSOl forms a disk-hke struc- 
ture, with a large amount of mass remaining in the high 
density inner region of the atmosphere. 

During the simulations, we trace the mass inflow 
through the inner boundary and the resulting total ac- 
creted mass is plotted against time in Fig. [7] as a fraction 
of the initial total cloud mass. Only the accreted mass 
originally in the cloud is taken into account. Shortly 
before G2's pericenter passage in 2013.5, material moves 
through the inner boundary for the first time, but it takes 
roughly until 2060 to swallow 40% of the cloud material 
in model CCOl. This corresponds to an additional av- 
erage accretion rate of 6 x 10^^ Mq yr"-'^ until the year 
2025. The time resolved accretion rate is shown in Fig.lH] 
for our standard cloud and ring model. During the whole 
simulation, the time averaged accretion rate remains ap- 
proximately constant. However, hydrodynamical insta- 
bilities break up the nozzle-like stream of gas leading to 
clumpy accretion. Hence we observe deviations from the 
mean of up to a factor of ten with a typical time scale 
of the order a few months. Concerning the ring model 
SSOl, the accretion through the inner boundary starts 
approximately half a year later and we estimate the av- 
erage accretion rate of cloud material to 4x 10~^ Mq yr~^ 
(until the year 2025), roughly 67 times higher compared 
to model CCOl, again remaining roughly constant on av- 
erage, but highly variable. However, these numbers have 
to be interpreted with great care: Firstly, they only con- 
cern cloud material and secondly, in these idealized sim- 
ulations, we do not take the back reaction of the cloud 
onto the ambient atmosphere into account. The latter 
could cause either enhanced accretion due to a fractional 
destabilization or a lower accretion rate due to the cre- 
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Figure 8. Measured mass accretion rate through the inner bound- 
ary of our standard model CCOl (black line) and the spherical ring 
model SSOl (red Une). 

ation of an even hotter inner bubble or feedback from the 
central source. What actually happens depends crucially 
on the detailed thermodynamics and the state of the hot 
atmosphere in this region, which is both not very well 
understood and not easily accessible via observations. 

5.4. Connection to observations and predictions for 
future observations 

Fig. [ol displays the position- velocity (PV) diagram of 
the whole orbit of G2. It shows the projected distance 
of the cloud material to the black hole on sky and the 
line-of-sight velocity of the cloud, taking the orbital an- 
gles as inferred from observations into account (Tab. Rj). 
Overlayed as filled contours are the data of the hydrody- 
namical simulation CCOl at various stages of its evolu- 
tion, as indicated in the plot. It represents the square of a 
mass histogram of the 2D gas distribution, representative 
for the Br 7 emission of the cloud gas. The cloud starts 
with overall constant velocity. When moving towards the 
SMBH, it develops a velocity shear, which is the largest 
during pericenter passage. After being disrupted shortly 




Figure 7. Accreted cloud mass of the models discussed in this 
paper relative to their initial cloud masses. See Table[2] for more 
details on the simulations. 



Figure 9. Position-velocity diagram for the simulation CCOl. 
The distance to Sgr A* — projected on the sky - is plotted against 
the line-of-sight velocity. The colored contours display the square 
of a mass histogram, relative to the total initial mass of the cloud, 
representative for the Br 7 emission. The dashed line represents 
the evolutionary path of a test particle in the center of the cloud 
on the orbit of G2 with the parameters given in Tablel!] for a full 
orbital revolution. The dotted lines show the result of the test 
particle simulation. 
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after the closest approach, the shear decreases when the 
cloud gets compressed again, but never reaches the val- 
ues expected from test particle simulations, as in these 
stages, hydrodynamical effects dominate the evolution 
of the cloud. Nonetheless, the maximum signal is in this 
model always expected to arise from close to the test par- 
ticle orbit of the cloud for the whole evolutionary time 

probed with our simulations. 

We zoom into the lower part of the diagram in Fig. [lO] 
for the CC model (panels a,b) and the SS model (pan- 
els c,d) in order to compare to available observations and 
make predictions for the near-future evolution. The com- 
pact cloud model was designed in order to give a good 
match with the observed head emission in the 2008.5 and 
2011.5 data, which is very well fulfilled. In contrast to 
this, we interpret the observations for the case of the 
spherical shell model such that the head (G2) and tail 
(G2t) emission are the focal points of an already dis- 
rupted large cloud complex (eG2). In order to get a 
reasonable match with the observed PV diagram, this 
requires a spherical shell of gas, realized as a ring in our 
two-dimensional simulations. As the ring starts at its 
apocenter distance of 0.05 pc at the inner edge of the 
disk(s) of young stars, it seems obvious to interpret it as 
being the resul t of stellar evolutionary processes. As can 
be seen in Fig 10 (panels c,d), it is in reasonable agree- 



ment with the combined observations of the head and 
tail emission of the cloud and can well explain the off- 
set of the tail from the inferred test particle orbit of G2. 
Fig. [Tl] show our predictions for observations in 2013. 
Distinct and observable differences are visible: Whereas 
our compact cloud simulation (panels a,b) produces emis- 
sion spread around the expected location of the test par- 
ticle orbit of the cloud, the spherical shell simulation dif- 
fers significantly by extending to higher velocities up to 
-1-2000 km s^i (panels c,d). 

In Fig. [12] we restrict the PV-diagram to the expected 
emission from the hydrodynamically formed tail of the 
cloud in the Compact Cloud Scenario. As expected, the 
stripped-off gas fills the interior of the test particle curve 
centered on G2 (thick black dashed line). The earlier 
the gas was removed from the compact cloud, the lower 
its present velocity due to the ram-pressure interaction 
and mixing with the ambient hot medium. The turbu- 
lent motions of this low-density material can be seen in 
the large width of the distribution in velocity direction, 
which even extends to velocities with the opposite sign. 
The contours vertically offset from the actual location 
of G2 could be related to the stripped-off gas from the 
front part of the cloud. Fig. 12 shows that what was 



often referred to as the tail oiG2 in the observed PV 
diagram, cannot be explained by stripped material from 
the compact cloud, causing problems for the Compact 
Cloud model or any model where the cloud is the result 
of a compact invisible so urce of gas ( Compact Source Sce- 
nario, see Sectionfl] and Murray-Clay & Loeb 2011). In- 
stead, the so far undiscussed contours at low distances to 
Sgr A* and low velocities could be related to a hydrody- 
namically stripped tail of gas. In contrast, our Spherical 
Shell Scenario is qualitatively able to explain the whole 
structure of the observed PV-diagram. 

A second promising possibility to distinguish be- 
tween different model pred ictions is via X-ray and IR- 
observations. In Sect. 5.3 we discussed the mass ac- 



cretion rate through t he inner boundary of our domain. 
Baganoff et al. ( 2003 ) observationally determine an ac- 
cretion rate of 10"*^ Mq vr~ ^ at the Bondi radiu s, located 
at 10^ Schwarzschild radii (Genzel et al. 2010). We use 
the scaling relation 



Afi, 



M, 



Bondi 



?^Bondi 



(7) 



where Hn, Bondi is the inner radius or the Bondi radius 
and Mjn^Bondi IS the mass accretion rate at the inner 
boundary o f our domain or the B ondi radi us, respec- 
tively (e. g. Blandford fc BegelmanlflQ QQ Igumenshchev 
et al.||2003[ ) and s = 0.27 ( | Yuan et arj|2003 ) to calculate 
the accretion rate expected at the inner boundary of our 
domain, which is located at roughly 1500 Schwarzschild 
radii distance from the SMBH. It results in a value of 
3 X 10~^ Mfr,yr~^. As discussed in Sect. 5.3 the mat- 



ter streams through our inner boundary in form of small 
droplets or filaments. Given the proximity of the SMBH, 
the evaporation timescale is very short (see equation p|. 
Hence our first assumption is that the cloud material will 
evaporate on its further way in and join the hot accretion 
flow. To estimate the effect on the observable luminosity 
of Sgr A*, we use the spectral energy distr i bution s de- 



rived for the ADAF models of Yuan et al. (2004). We 



interpolate their SEDs for various mass accretion rates 
shown in their Fig. 5 to match the average mass accretion 
rates of the flow through the inner radii of our models 
and compute the IR (0.8 — 2.5/im) and X-ray (2-8 keV) 
luminosity. If the mass accretion proceeds in this mode, 
we do not expect a significant boost of the X-ray or IR 
luminosity for the CCOl model, as even the spikes ten 
times above the average accretion rate would only lead 
to a factor of 1-2 increase of the IR and X-ray signal. 
In contrast to this, model SSOl would result in an X-ray 
luminosity of 4 x lO'^^ergs^^, corresponding to a boost 
of a factor of roughly 80 compared to the quiescent X-rav 



luminosity of Sgr A* (4 x 10 ergs derived from Yuan 



et al. 2004 



to 5 X 10' 



Fig. 5). The expected IR luminosity amounts 
'ergs~^, roughly the same factor above the 
assumed quiescent value (7 x 10"^^ erg s^^). Concerning 
the spikes of the accretion rate distribution (roughly a 
factor of ten above the average accretion rate), we cal- 
culate boost factors of roughly 3 x 10* in the X-ray and 
6 X 10'^ in the IR. These high boost factors would enable 
us to test the validity of our spherical shell scenario in 
the coming years. However, another possibility (we are 
unable to test with our current models) is that the cold, 
disrupted cloud significantly affects the thermodynamic 
structure of the hot accretion flow close to the GC. This 
might then prevent the complete thermal evaporation of 
the cloud fragments and enable the formation of a thin 
and relatively cold accretion disk. To get a first rough 
estimate of the expected signal, we use the average mass 
transfer rate through the inner boundary of our domain 
(Fig. Is]) and assume it is representative for the steady- 
state mass accretion rate through such an inner accre- 
tion disk, which we assume forms instantly. Given this 
assumption, the bolome tric lu minosity of the disk can be 
estimated to be (Krolik..l999): 
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Figure 10. Comparison of the position- velocity diagrams for the simulations CCOl (panels a, b) and SSOl (panels c,d) with observations 
(background contours). The distance to SgrA* - projected on the sky - is plotted against the line-of-sight velocity. The colored contours 
display the square of a mass histogram, relative to the total initial mass of the respective cloud model. The dashed line represents the 
evolutionary path of a test particle in the center of the cloud on the orbit of G2 with the parameters given in TablelT] for a full orbital 
revolution. The green dotted lines show the result of the test particle simulation. 
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Figure 11. Predictions for the distribution of the compact cloud model CCOl (panels a, b) and the spherical shell model SSOl (panels c,d) 
in the position- velocity diagram for the year 2013. The distance to SgrA* - projected on the sky - is plotted against the line-of-sight 
velocity. The colored contours display the square of a mass histogram, relative to the total initial mass of the respective cloud model. The 
dashed line represents the evolutionary path of a test particle in the center of the cloud on the orbit of G2 with the parameters given in 
Table[T]for a full orbital revolution. 

lication, but might influence these values significantly. 
However, a similar increase of the luminosity compared 
to the curren t state might in tu rn influence the RIAF so- 
lution itself (Yuan et al. 2004) as well as the dynamical 
evolution of the cloud via its radiation pres sure interac- 
tion, as has fo r example been investigated by |Schartmann 
et al. (2011). More detailed modelling of the physical 
processes at work as well as the hot atmosphere is nec- 
essary to determine the expected light curve, triggered 
by the future evolution of the cloud. Thermal conduc- 
tion might be one of the key players here as the cloud 
evaporation timescale becomes very short in this phase. 



This results in an average bolometric luminosity of 
8x10^ Lq and 5x10^ Lq for the cloud model and the ring 
model respectively (averaged over a time period from the 
start of the accretion until the year 2025). In order to 
derive the observable X-ray luminosity, we use the ac- 
cretion disk spect r al ene rgy distribution constructed in 
Schartmann et al. ( 2005 ) from a combination of observa- 
tions and simulations. I'his results in averaged values of 
i2-8kev = 4x lO^^ergs"^ and L2-8keV = 3x lO^^ergs"^ 
for the cloud and ring model in the energy range between 
two and eight keV and 3 x lO'^^erg s~^ and 10'^^ erg s^^ for 
the two models in the IR, the latter being orders of mag- 
nitude above the current quiescent luminosity. Which 
fraction of the cloud material will accrete in this mode 
is currently unclear and the detailed prediction of the 
resulting SED is clearly beyond the scope of this pub- 



5.5. Influence of the equation of state 

A test simulation has been run with an adiabatic equa- 
tion of state (simulation CC04). This represents ineffi- 
cient external cooling and heating sources, but takes adi- 
abatic cooling and heating via expansion and contraction 
into account. The early evolution of the density distri- 
bution is very similar in both cases. Major differences 
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0.4 0.3 0.2 0.1 

offset from SgrA« [ob] 

Figure 12. Position-velocity diagram for the simulation CCOl, 
where we only show the emission of the tail of the cloud, selected 
by tracer values between the threshold value and 0.5. The distance 
to Sgr A* - projected on the sky - is plotted against the line-of- 
sight velocity. The colored contours display the square of a mass 
histogram, relative to the total initial mass of the cloud. The thick 
dashed line represents the evolutionary path of a test particle in 
the center of the cloud on the orbit of G2 with the parameters 
given in TablelTlfor a full orbital revolution. 

only arise after the pericenter passage. A comp aris on of 
the state in the year 2030 is presented in Fig. [13] We 
find that in simulation CCOl, the pressure inside the re- 
maining cloud gas is larger compared to the adiabatic 
case CC04. The reason is that the gas in the adiabatic 




simulation is able to cool when the cloud expands and 
gets disrupted after the pericenter passage. This effect is 
also visible in the formation of a low density but high 
pressure region in downstream direction of the cloud, 
which is less pronounced in the adiabatic simulation. The 
lower temperature and pressure in the adiabatic case lead 
to a more filamentary cloud with a stronger density en- 
hancement, whereas in the other case, the larger pressure 
support leads to a more uniform density distribution, 
spreading over a larger area. This example shows again 
that the early evolution of the cloud can easily be under- 
stood with simple hydrodynamical simulations and we 
can trust the predictions for the near-future observations. 
However, the late-time evolution - after tidal disruption 
has severely affected the cloud structure - significantly 
depends on the detailed physical processes and the struc- 
ture of the surrounding atmosphere. For an in-depth 
understanding of the long-term fate of the cloud, more 
sophisticated (magneto)-hydrodynamical simulations are 
required. 

5.6. Influence of the atmospheric density distribution 

Our knowledge of the temperature and density distri- 
bution of the hot atmosphere in the Galactic Center re- 
gion is based on theoretical models, which are consistent 
with the small number of data points observed with the 
Chandra satellite and the rotation measures (see above). 
Hence, it is interesting to ask whether the current and fu- 
ture evolution of the cloud will enable us to probe the hot 
atmosphere. A first impression of the effect of the den- 
sity distribution of the ambie nt at mosphere on the cloud 
evolution was given in Sect. 



2.2 As was shown there 



Figure 13. Comparison of the density of the cloud in 2030 for our 
standard model (panel a, CCOl) and the adiabatic test simulation 
(panel b, CC04). 



under the assumption of initial pressure equilibrium, the 
starting date would be changed by five to ten years if the 
atmospheric density distribution would be changed by a 
factor /hot of 0.5 (CC05) or 2.0 (CC06). However, the ab- 
solute values of the preferred origin date of the cloud de- 
pend on the not very well known total mass of the cloud 
and hence our results should be taken with care. Here, we 
briefiy discuss the hydrodynamical realizations of these 
two best-fit models from the test particle simulations for 
the Compact Cloud Scenario. As expected, both clouds 
are in good agreement with the observed PV diagrams 
in 2008.5 and 2011.5, but the late-time evolution of the 
density distribution changes significantly, as can be seen 
in Fig. |14| The growth of hydrodynamical instabilities 
scales wrth the density contrast of the cloud/atmosphere 
system and the importance of ram pressure effects in- 
creases with the density of the environment (see Sect.B. 
The earlier onset of KHI in the high-density case -al- 
ready during the pericenter passage - leads to a more 
filamentary cloud boundary. The interior is affected as 
well due to the partly disruption and the later buckling 
of the cloud. Gas transfer through the inner boundary 
starts roughly at the same time in both cases. However, 
the reduced ram pressure for the case of simulation CC05 
leads to a smaller accretion rate (blue graph in Fig. JTl 
compared to the high density atmosphere (CC06), which 
displays a steeper increase of the accretion rate with time 
(yellow curve in Fig. I7| . Concerning the evolution in the 
PV-diagram, model T!C06 spreads over a larger region 
during the late evolution due to the faster disruption of 
the cloud in the high density case. 
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Figure 14. Comparison of two compact cloud models with a factor of two decreased density of the atmosphere (panel a,b, CC05) and a 
factor of two increased density of the atmosphere (panel c,d, CC06) with respect to our standard model. 



6. DISCUSSION 

Numerically, the most critical point of the simulations 
is the treatment of the atmosphere. As already discussed 
in Sect.[3J the ADAF solution which fits the Chandra ob- 
servations of the hot atmosphere best is unstable to con- 
vection and cannot be modeled in a simple way. In order 
to enable a first investigation of the origin and further 
evolution of the cloud, we artificially stabilize the at- 
mosphere and make several simplifying assumptions: (i) 
We neglect the accretion flow of the diffuse atmosphere 
towards the center and set up the central hot gas distri- 
bution in hydrostatic equilibrium, (ii) We suppress the 
growth of disturbances in the atmosphere by re-setting 
those cells within the atmosphere, which have not inter- 
acted with the cloud yet. This is necessary in particular 
close to the central boundary condition, where distur- 
bances due to the Cartesian grid start to increase due to 
the entropy gradient. To avoid these complications, we 
use a passive tracer field. Cloud material is lost in those 
regions of the computational grid where the fraction of 
cloud gas relative to atmospheric gas is below a thresh- 
old of 10""'. This leads to a roughly linear increase of 
the mass loss starting around the time of the pericenter 
passage in the standard model CCOl. At the end of the 
simulation in the year 2060, roughly 1.3% of the initial 
cloud mass is removed from the simulation. This also 



means that the angular momentum and kinetic energy, 
which the cloud transfers to the atmosphere is finally 
taken out of the simulation. A critical discussion of the 
infiuenc e of the threshold value for the tracer is given 
in Sect. 16.21 Further limitations of our simulations are 
the negligence of thermal conduction, which will be im- 
portant at l east in the late phase of the evolution of the 
cloud (Sect. 4.3) and the restriction of the simulations to 
two dimensions. Perpendicular to the modelled orbital 
plane, we expect the cloud to collapse under the influ- 
ence of the gravitational force of the black hole, forming 
a thin disk and later a thin filament close to pericenter. 
After passing the black hole, the cloud is expected to 
expand again in vertical direction due to the overpres- 
sure with respect to the surrounding medium. Civen the 
larger contact surface with the hot ambient medium in 
rest, enhanced stripping due to KHI is expected, which is 
most important for the late time shredding of the cloud 
after pericenter passage. In contrast to this, we expect 
less ram pressure interaction in the 3D case, which might 
affect the angular momentum exchange with the atmo- 
sphere and hence slightly alter the mass transfer towards 
the center. 
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Figure 15. Resolution study of our standard model. Displayed 
is our standard resolution (panel b, CCOl), half the resolution 
(panel a, CC07) and twice the resolution (panel c, CC08). 

6.1. Resolution effects 

The effects of resolution on the density evolution of 
our simulations are displayed in Fig. |15| The snapshots 
are taken from the compact cloud simulation starting on 
G2's orbit in the year 1995.5. Panel b displays our stan- 
dard resolution, whereas we decreased the resolution in 
both axis-directions by a factor of two in panel a and 
increased the resolution by a factor of two for the sim- 
ulation displayed in panel c. The snapshots correspond 
to the year 2016, shortly after the pericenter passage. 
The basic dynamical evolution as well as the disruption 
due to tidal forces is very comparable. However, the ex- 
pected resolution effects are visible in the detailed evo- 
lution of the hydrodynamical instabilities. This leads to 



a stronger disruption of the nozzle-shaped inflow of gas 
towards SgrA* with increased resolution. In the high- 
est resolution simulation, the stream has broken up al- 
ready and the mass-flux through the inner boundary is 
clumpy, compared to the smooth accretion of gas through 
the nozzle in our low resolution run. Our standard reso- 
lution run already shows some clumpy accretion, which 
also increases at later stages of the evolution (Fig. |4] 
lower row). The mentioned accretion rates through the 
inner boundary differ by about 10% between the highest 
and the lowest resolution in our study. 

6.2. Influence of the choice of the tracer threshold value 

Fig. [l6] displays the influence of the tracer thresh- 
old value. This is done for our simulation starting on 
G2's orbit in the year 1995.5 (model CCOl) and we 
analyse the state of the simulation in the year 2025.5, 
as only the late-time evolution shows significant differ- 
ences. From panel a to c, the tracer threshold value in- 
creases from 10~^ (panel a) to lO"** (panel b) and 10^^ 
(panel c). Again, as in the resolution study, the basic 
dynamical evolution is unaffected by the choice of the 
threshold value. However, in this series, two system- 
atic effects are visible: (i) For too high threshold values 
(Fig. 16:), the shape of the outer boundary in upstream 
direction changes slightly. As soon as the cloud mixes 
with the surrounding medium, gas is removed from the 
simulation. This artificially suppresses instabilities, (ii) 
For too low threshold values, the disturbed atmosphere 
in downstream direction of the cloud gets susceptible to 
the growth of convectively unstable regions. Bubbles are 
formed behind the cloud as well as in the disturbed atmo- 
sphere along the traversed orbit of the cloud. Quantify- 
ing the effects of the de-stabilization of a realistic atmo- 
sphere and its back-reaction onto the further evolution 
of the cloud is beyond the scope of these idealized simu- 
lations. Given the considerations briefly discussed here, 
we decided to use a tracer threshold value of 10~^. 

6.3. Possible influence of magnetic fields on the 
evolution of the cloud 

Strong magnetic fields of the order of 100 fiG (with a 
lower limit of 50 /iG on 400 pc scales) have been observa- 
tionally inferred near the Galactic Center JYu sef-Zadeh 
fc Morris|19"87l|Morris fc Yuse f-Zadch"1989| |Crocker et at 



20ldr Field strengths of this magnitude might signifi- 



cantly influence the dynamics and physical interactions 
occuring in the interaction region of the cloud. Several 
effects might change the dynamics and appearance of 
the cloud: (i) magnetic pressure inside the cloud and of 
the hot atmosphere changes the total pressure balance 
and hence the size along the cloud's orbit, (ii) magnetic 
tensi on forces can lead to a conflnement of the cloud 
(e. g. Krause et al. 2012 ), thereby alleviating the problem 
of the too fast tidal disruption, (iii) depending on field 
morphology, the Kelvin-Helmholtz-Instability can be sig- 
nificantly suppressed, again leading to a longer survival of 
the cloud and significant changes in the late time evolu- 
tion and (iv) magnetic fields lead to anisotropic thermal 
conduction, altering the cloud evaporation time scale. 
Hence it will be necessary to consider magnetic fields in 
future simulations determining the fate of the observed 
cloud. The problem is that this adds additional degrees 
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Figure 16. Studying tiie influence of the tliresliold value of the 
tracer field for artificially stabilizing the atmosphere. Displayed is 
a value of 10"'' (panel a), 10~^ (panel b, our standard model) and 
10-2 (panel c). 

of freedom as neither the strength of the fields nor its 
morphology is currently known in the sub-parsec sur- 
rounding of Sgr A*. 

7. CONCLUSIONS 

In this paper, we present a set of idealized hydrody- 
namical simulations to study the origin and the fate of 
the recently discovered cloud complex in the Galactic 
Center. Building up on analytical estimates in paper I, 
we concentrate on two possible scenarios: (i) the Com- 
pact Cloud Scenario, where the cloud is assumed to 
have formed in the recent past on the orbital path of 
the G2 component of the cloud and (ii) the Spherical 
Shell Scenario, in which the observations are inter- 
preted as an already disrupted spherical shell (or ring, 
respectively in two dimensions). The early evolution un- 
til close to the pericenter is dominated by tidal interac- 
tion in both cases. Shortly after pericenter passage, when 
the cloud is already tidally disrupted, secondary hydro- 
dynamical instabilities as well as ram pressure become 
important, leading to some angular momentum redistri- 
bution and accretion towards the center. Nonetheless, 
the cloud keeps a large fraction of its initial angular mo- 
mentum, which efficiently hinders direct accretion onto 



Sgr A*. Accretion is expected to occur via a thin nozzle- 
like feature connecting the cloud with the SMBH, as a 
result of continuous ram-pressure interaction of the up- 
stream boundary of the cloud with the dense inner part 
of the atmosphere. Hence ~ depending on the destabiliza- 
tion of this stream (e. g. by thermal conduction) - there 
might be enhanced cold gas accretion within the next 
years, which might boost the luminosity of Sgr A* in an 
observable way. However, according to our simulations, 
a large fraction of the cloud will remain on a similar orbit 
compared to the test particle orbit and lead to a steady 
infiow of gas towards the Galactic Center. In contrast to 
the early evolution of the cloud, the post-pericenter stage 
depends sensitively on the detailed modelling of physical 
processes as well as the atmospheric density distribution, 
necessitating more detailed hydrodynamical simulations 
to predict the (far) future evolution of the cloud com- 
plex. A more realistic treatment of the unstable nature 
of the atmosphere might influence our derived accretion 
rates substantially, when taking the back-reaction of the 
cloud into account. 

In summary, we currently favor the Spherical Shell 
Scenario for the following reasons: (i) It allows the cloud 
to be started within the range of the disk(s) of young 
stars, (ii) Being a spherical shell of gas further strength- 
ens the notion of it being the result of stellar evolutionary 
processes, (iii) The model is able to explain the major 
part of the structure (eG2) seen in the observed position- 
velocity diagram, whereas the Compact Cloud Scenario 
only captures the dominant part of it (G2). (iv) Tidally 
stripped gas cannot account for the G2t component of 
the observed PV-diagram, which disfavors the Compact 
Cloud Scenario and any Compact Source Scenario, but a 
so far undiscussed component of the PV-diagram could 
be directly related to it. 

We are grateful to Romain Teyssier for helpful discus- 
sions and an anonymous referee for helpful comments. 
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